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Abstract 

The dynamics of a qubit in a structured environment is investigated theoretically. One point of 
view of the model is the spin-boson model with a Lorentz shaped spectral density. An alternative 
view is a qubit coupled to harmonic oscillator (HO), which in turn coupled to a Ohmic environ- 
ment. Two different methods are applied and compared for this problem. One is a perturbation 
method based on a unitary transformation. Since the transformed hamiltonian is of rotating wave 
approximation (RWA) form, we call it the transformed rotating wave approximation (TRWA) 
method. And the other one is the numerically exact method of the quasi-adiabatic propagator 
path-integral (QUAPI) method. TRWA method can be applied from the first point of view. And 
the QUAPI method can applied from both points of views. We find that from the 1st point of 
view QUAPI only works well for large T. Since the memory time is too long for the practical 
evaluation of QUAPI when T is small. We call this treatment as QUAPI1. And from the 2nd point 
of view, QUAPI works well for small T, since the non-adiabatic effect become more important as V 
increases, one need smaller time-step and more steps to obtain accurate result which also quickly 
runs out the computational resources. This treatment is called QUAPI2. We find that the TRWA 
method works well for the whole parameter range of T and show good agreement with QUAPI1 
and QUAPI2. On the other hand, we find that the decoherence of the qubit can be reduced with 
increasing coupling between HO and bath. This result may be relevant to the design of quantum 
computer. 
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I. INTRODUCTION 



Dissipative quantum dynamics is of crucial interest among scientist. Since the quantum 
dynamics are always inevitably affected by its environment, various physical and chemical 
phenomena are related to the dissipation, range from the spontaneous emission to electron 
transfer in molecular, from qubit decoherence to photon harvest in photosynthesis l|-|5]. 
Spin-boson model, the simplest possible model to describe dissipation, offers a comprehen- 
sive way to study the decoherence phenomenon. In the context of the electron transfer in 
molecular, spin boson model has been studied intensively over the past decades, and it shows 
revival interest among scientists because of the possible application of quantum computation 
and information. 

For the spin boson model, the environmental property is characterized by the spectral 
density J(u), which is usually assumed to be a power law distribution, J(w) oc u s , it is 
called sub-Ohmic when < s < 1, Ohmic when s = 1 and super-Ohmic when s > 1. The 
most studied case is the Ohmic spectral density, it describes a case when the dissipation 
is the same at all frequencies, which is the case for the many environments. A physical 
example corresponds to the Ohmic case is the dissipation in a pure resistor circuit. For the 
sub-Ohmic case, it arouse lots of interest recently, because of the controversy related to the 
quantum phase transition. 

In this work, we study the dynamics of a qubit in a structured environment. Two points 
views are available for the problem we are interested in. One point of view is the spin 
boson model with a Lorentz shaped spectral density. And the alternative view is a two- 
level system (TLS) coupled to harmonic oscillator (HO), which in turn coupled to a Ohmic 
environment. The model describes some situations in experiments. For example, a flux qubit 
is usually read out by a dc-SQUID with a characteristic plasma frequency, consequently, the 
environmental noise of the SQUID is also transferred to the qubit leading to decoherence 
and dissipation. It also describes a Cooper-pare box (CPB) coupled to a transmission line 
resonator or a qubit placed in a leaky cavity. 

From the first point of view, many tradition treatment of spin-boson model can be applied, 
some of them are largely numerical, such as the quantum Monte Carlo 



17] . real-time 



JL5 

renormalization group [18| : _ quasi-adiabatic path integrals (QUAPI) |19|-22|. low equation 



renormalization (FER) 
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28j and numerical renormalization group (NRG) 
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j_others are mainly analytical such as the non-interacting blip approximation (NIBA) 
434j. rigorous Born approximation 25|,|35| or Bloch-Redfield 
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371 ] , and the transformed 



rotating- wave approximation (TRWA) method 38N41L Ti 



density, it has been studied by FER 



26|-|28|, NIBA"j27,|33|, Bloch-Redfield 



1 now, for this particular spectral 



27J and TRWA 

40l . |41| . And from the second point of view, it has been studied by QUAPI (22) , the van Vleck 



perturbation theory together with a Born-Markov master equation (VVBM) 



equation based on a single mode TRWA 



42|, Redfield 



43| . Although it has been studied intensively, the 



most studied case is the weak HO-bath coupling. While in this work, we explore the qubit 
decoherence in both the weak and strong HO-bath coupling regimes. 

One commonly believed concept is that temperature and the coupling to noisy bath only 
olay a negative role in preserving the qubit coherence. However, it is pointed out in Ref. 



44J that the temperature can help the coherence when the qubit is coupled to a TLF (or 



spin-boson) environment. We would like to ask: Is it possible to reduce the decoherence by 
increasing the system-bath coupling? In this paper, we examine the coupling effect on the 
quantum decoherence, we study aforementioned problem from the 1st point of view by the 
TRWA method in both the weak and strong HO-bath coupling regimes. To check the result, 
we also use the numerically exact method of the quasi-adiabatic propagator path-integral 
(QUAPI) method. We find that from the 1st point of view QUAPI only works well for large 
T. Since the memory time is too long for the practical evaluation of QUAPI when T is small. 
We call this treatment as QUAPI1. And from the 2nd point of view, QUAPI works well 
for small T, since the non-adiabatic effect of the Ohmic bath become more important as V 
increases, consequently, smaller time-step is required, thus, for the same memory time, one 
need more steps to obtain accurate result which also quickly runs out the computational 
resources. This treatment is called QUAPI2. We find that the TRWA method works well 
for the whole parameter range of T and show good agreement with QUAPI 1 and QUAPI2. 
On the other hand, we find that the decoherence of the qubit can be reduced with increasing 
HO-bath coupling when the HO-bath coupling is larger than the qubit-HO coupling. 

The paper is organized as follows. In sec. UH we present explicitly the two points of views 
of the problem we are interested in. In sec. II I II and IIVI we briefly introduce the TRWA 
and the QUAPI method, and discuss how to adopt these methods. In sec. |V] the coupling 
effect on the quantum coherence is discussed and a comparison of the TRWA with QUAPI 
is given. Finally, a brief conclusion is presented in Sec. IVII 
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II. TWO DIFFERENT POINTS OF VIEWS OF THE PROBLEM 



The 1st point of view of the problem is the spin-boson model (un-biased case) 

H = -yas +J2 u kblh + \a*Y^9kQ>k + b k), (!) 

k k 

where a x and a z are the pauli matrices, A is the gap of of the qubit, bk (or b\) are the 
annihilation (or creation) operators of the bath. And the spectral density is Lorentzian, 

J(U) = E ^ - Uk) = /02 ,,^lr..n^ ( 2 ) 



(n 2 



This spectral density possesses a shape peak at position Q especially when T is small, which 
may challenge the conventional method of the spin boson model. In the limit of Q — > oo, 
J{oj) is reduced to the Ohmic spectral density. 

The above model can be exactly mapped to a model where the qubit is dissipated by a 
multi-mode Ohmic bath via a harmonic oscillator, which is the 2nd point of view of this 



problem 



45 



46| . The Hamiltonian reads (H — 1) 



+ (B^ + B) 



9o a z + ftfc(ftfc + h 



k Uk 



(3) 



where, B (or B^) are the annihilation (or creation) operators of the HO, bk (or b\) are the 
annihilation (or creation) operators of the corresponding bath mode. The last term is the 
counter-term, which cancels the additional contribution due to the coupling of the HO to 
the bath 0, l^] . The corresponding spectral density is of Ohmic form 

J(u) = K l 6 ( u - S fe) = To;e- w/aJc (4) 
k 

where u c is the cut-off frequency (throughout this work, we use the value u c = 20ft). The 
relation between g and a follows as g = ^v/gf ■ From this 2nd point of view, the over all 
system TL plus HO can be considered as the system which is then dissipated by an Ohmic 
bath. 

III. TRANSFORMED ROTATING- WAVE APPROXIMATION METHOD 



The transformed rotating-wave approximation 
and then applied to a sequence of problems 
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method is developed by one of the author 
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4601 ] Our starting point is the SBM 



with the structured spectral density. A unitary transformation is first applied to the SBM 
Hamiltonian, H' = exp(S)H exp(— S), with the generator S = J2k ~ bk)cr z . The 

transformed Hamiltonian can be decomposed into three parts: 

fc k " 



H[ = ^J2^( 1 -^M + h)--fr j AX, (6) 

H' 2 = -yA {coshX- 77} - ^AjsinhX- r?xj, (7) 

where, X = J^ fe ^^k{b\ — and is the thermodynamic average of coshX. In the limit 
of zero temperature it is, 



rj = exp 



a 2 



2u 



2^k 
k 



(8) 



Obviously, H' Q can be solved exactly since the spin and bosons are decoupled in H' . r/A gives 
a rough approximation of the renormalized qubit frequency and (rj — 1) A is the corresponding 
Lamb shift of the qubit due to the coupling to the bath. The eigenstate of H' can be 
expressed as direct product: |s)|{nfc}), where \s) is the eigenstate of a x and ({n^}) is the 
eigenstate of phonons, which means that there are n k phonons for mode k. Therefore, the 
ground state of H' Q is given by : \go) = | si)|{0fc}), where | si) is the lower eigenstate of spin 
and |{0jfc}) stands for the vacuum state of the bosons. 

The choice of rj in Eq. ([8]) insures H' 2 contains only the terms of two-boson and multi- 
boson non-diagonal transitions and its contribution to physical quantities is (gl) 2 and higher. 
Therefore, H' 2 can be omitted in the following discussion. If we let H[\go) = 0, which 
ensures H[ to be small and more suitable for the subsequent perturbation treatment, then 
the parameters £fc's are determined as, 

U k +T]A 

Consequently, H[ is transformed to the rotating-wave form, 

k 

where V k = rjAg k ^ k /uj k = g k rjA/{u k + t]A) and a± = (a z ^fia y )/2. Actually, some effect of 
the anti-rotating-wave terms has been taken into account in H[ by the unitary transformation 



which is embodied in the renormalized coupling constant Note that ~ 1 if the boson 
frequency u k is larger than the renormalized tunneling r)A, but <C 1 for <C r?A. Since 
the transformation generated by S is a displacement of bosons, physically, one can see that 
high-frequency bosons (uk > i]A) follow the tunneling particle adiabatically because the 
displacement is gk^k/^k ~ 9k/^k- However, bosons of low-frequency modes u>k < i]A in 
general are not always in equilibrium with the tunneling particle, hence the particle moves 
in a retarded potential arising from the low-frequency modes. When the non-adiabatic effect 
dominates, Uk <C r]A, the displacement «C 1. 

The total density matrix (system+environment) x'{t) obeys the Liouville-von-Neumann 
equation, 

j t X'(t) = -z[H[(t), X '(t)}, (10) 

where the tildes denote operators in the interaction picture with respect to H' . Iterating up 
to the second order and tracing out the environmental degrees, one get the master equation 
within the Born approximation, 

j t p'(t) = -iTQ: s [£ri(t), <S> ^(0)] - dt^ B [^(t), [^Ti(0, A ® ^(0]], (11) 

where p' is the reduced density matrix p'{t) =Tr B [%'(£)], and x'(t) is replaced by an ap- 
proximate factorized density matrix x'(t) ~ p'b® P'(t)- The environment is usually assumed 
to remain in thermal equilibrium p' B = e~^ H ' B /Tre _/3H s with H' B = ^2cUkb k bk, which is jus- 

k 

tified when the environment is 'very large' and the coupling 'weak' (Vk <C A, Q), so that the 
back-action of the system onto the environment can be neglected. Substitute H[ into the 
master equation, we get 

^ = -iTr B [H[,m - $> fc 2 fdt'XM (12) 

k ^° 

where, 

X(M') = n fc [«r_(t)o- + (0/5'(0-o-+(Op'(0^-(*)]e <w * (t - t ' ) 

+ (n k + l)[e + (t)e_{t>)p'{t>) - a_{1f)?(ir)a + {t)]e-**W 
+ (n fc + lMi>+(*>-(0 " ^{t)~pf{t')a + {t')\e^-^ 

The master equation Eq. (Tl2l which is a 2 x 2 matrix equation can be solved exactly by 
the Laplace transform since the convolution theorem can be applied to the equation of each 



matrix element. Here, for simplicity, we only present the comparatively brief expression. 
Suppose the system is in the upper eigenstate of a z at time t — 0. At the zero temperature, 
the population difference P(t) = (<J z (t)) = Tr A (a z p(t)) is evaluated as 

p { t)= i -r d u iM^p. . (13) 

Wo [u - V A - R{u)} 2 + 7 2 (w) 
R(u) and 7(0;) in Eq. ( fl3l) are the real and imaginary parts of J2k 7 ( u ~ ^ + ~ u k) 

7 (w) = irJ(uj)(riA) 2 /(u + riA) 2 . (15) 



IV. QUASI- ADIABATIC PATH-INTEGRAL METHOD 



The QUAPI method is a numerical scheme based on a exact methodology 19M22I]. The 
starting point of QUAPI method is the generic system-bath Hamiltonian 

P 2 1 

j 3 

where, Hq is the Hamiltonian for the bare system, s is the system coordinate, and Qj 
are harmonic bath coordinates which are linearly coupled to the system coordinate. The 
characteristics of the bath are captured in the spectral density function 

j( W ) = |£_i-tf( W - Wi ). (i6) 

The reduced density matrix of the system evolve as p(s", s', t) = 
TTbath (s"\ e~ tHot p(0)e lHot \s'). If the path integral representation is discretized by N time 
steps of length At = t/N and the initial density matrix is assumed to be p(0) = p s (0)pbath(^), 
the reduced density matrix takes the form 

P (s", ^) = EE-EEEE ^ e ~ lHoAt • • • Wl e ' iHoAt \ s o) 

S N-1 S N-1 S l S l S S 

( S o \ Ps(0) | s ) ( S I e ° | S 1 ) ' ■ ■ ( S iV-l| e ° I s ) l( s 0i s ) s l~) s l)"' > S V-1> S 7V-1> S ! S ' ^0 

(17) 

where the discrete variable representation (DVR) is used, the symbol (k = N — 1) 

denotes the system coordinate at the time kAt on the forward and backward discretized 
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Feynman path. \s k } (k = N — 1) are the eigenstates of the system coordinate op- 
erator s. If a symmetric splitting of the time-evolution operator is employed e~ lHAt = 
e -tH env At/2 e -iH At e tH env At/2 w ^j 1 u ena = if — H , the corresponding influence functional 

reads 



Tr, 



bath 



p iH env (s")At/2 e -iH env {s+^ 1 )At _ _ _ e -iH env {s+)At/2 



Xp bath (0)e iHenv ^ a S > At / 2 ■ ■ ■ e iH env{s N _ 1 )At e iH,, nv {s')At/2 ^ ^g>j 

One can find that the equilibrium position of the bath mode is adiabatically displaced 
along the system coordinate. If Hq provides a reasonable zeroth-order approximation to the 
dynamics, the quasi-adiabatic propagator is accurate for fairly large time steps. That is 
the quasi-adiabatic partitioning is a good representation when the bath property is mainly 
adiabatic, where the bath can keep up with the motion the system quickly. And the discrete 
path is to take into account of the non-adiabatic effect. For most of case, the quasi-adiabatic 
partitioning is reasonable especially when the system bath coupling is not strong. Therefore, 
the QUAPI discretization permits fairly large time steps when the adiabatic bath dominates 
the system dynamics. If the bath is purely adiabatic, even no discretization is needed. In the 
continuous limit (that is for At — > 0, iV — > oo) the influence functional has been calculated 
by Feynman and Vernon 

l~J*dlf J* dt" [s + {t') - 8 -(t) 
x [a(t' - t")s + (t") - a* {t' - t")s-(t")] 

3 J ) 

where a(t) is the bath response function, which can be expressed in terms of the spectral 
density as 

1 f°° 
a(t) = — duJ(u) 
K Jo 



exp 



, . fiujh , . . . , . 

coth I — - — cosyojjt) — i sm(ujjtj 



(20) 



The last term in Eq. (ll9j) arises from the "counter-terms" which are grouped with the bath 
Hamiltonian in the quasi-adiabatic splitting of the propagator. With the quasi-adiabatic 
discretization of the path integral, the influence functional, Eq. [191 takes the form 

f 1 N k 

I = exp I -- Yl i s t ~ s k] [Vkk'st, - r)lk> s k>] 

{ k=0 fc'=0 



where = s" and sjf = s'. The coefficients r] kk / can be obtained by substituting the 
discretized path into the Feynman- Vernon expression Eq. (1191) . which is given in Ref. 201 ] . 

The QUAPI method is essentially a tensor multiplication scheme, which exploits the 
observation that for environments characterized by broad spectra the response function 
a(t) decays within a finite time interval. From the expression of the Feynman and Vernon 
influence funcitonal Eq. (FT?]) , one can see that a(t) characterizes nonlocal interactions, 
which connects system coordinate s(f) with s(t"). The path s ± (t') at time t' is connected 
to the all the paths s~(t") at earlier times, which makes the evaluation of Eq. (ITT)) a hard 
task. However, for a bath with a broad spectral density, such as a power law distribution of 
the spectral density, a(t) has the finite memory, the memory length typically extending over 
only a few time slices when the quasi-adiabatic propagator is used to discretize the path 
integral. After discarding the negligible "long-distance interaction" with t' — t" > Ak max At 
(or k — k'> AA; max ), the resulting path integral can be evaluated iteratively by multiplication 
of a tensor of rank 2Ak max . In other words, there exists an augmented reduced density tensor 
of rank Ak max that obeys Markovian dynamics. The details of the multiplication scheme is 
discussed to a great extent in the literature, here we only present the essential parameters 



and mention briefly how to adopt it to our specific problem 
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422|]. 



A. QUAPI1: treatment from the 1st point of view 

The Quapi method is based on an exact methodology, and for most conventional bath, 
such as bath with power law distributed spectral density, QUAPI work efficiently and con- 
verge easily. However, for the spectral density studied here, it pose challenge to the ap- 
plication of the QUAPI method, since the memory length is too long for the practical 
implementation of the QUAPI when F is small. The response function is depicted in Fig. 1, 
one can find that because of the characteristic peak frequency Q of the spectral density, the 
response function possess a coherent oscillation with frequency Q, and possesses very long 
memory time. However, when T becomes large, the damping of the kernel becomes signifi- 
cant, and the response function quickly damp to the zero, which enables the application of 
the QUAPI. Here we discuss briefly the parameters used in the QUAPI1 method: 

(i)The first parameter time-step At used for the quasi-adiabatic splitting of the path- 
integral. The memory time of the non-Markovian steps used by QUAPI is Ak max At. The 
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stability of the iterative density matrix propagation ensures the choices of At, it should 
not be too big nor too small, since the non-adiabatic effect requires more splitting of the 
path integral, that is smaller At. Whereas, since the memory length Ak max At is usually 
a fixed value for a particular bath, QUAPI method prefers larger At, and consequently 
smaller Ak max in consideration of the numerical efficiency (note that the algorithm scales 
exponentially with A/c max , also see the discussion of the second parameter Ak m3X ). Therefore, 
we should choose appropriate At to take into account both the non-adiabatic effect which 
prefer smaller time splitting and the non-Markov effect which prefer long memory time, 
typically, we choose At around J^-, that is to choose tens of fraction of the cycle time of 
the bare system dynamics. 

(ii) The second parameter is the memory steps Ak max . If Ak max < 1, the dynamics is 
purely Markovian. If the non-locality extends over longer time, terms with Ak max > 1 have 
to be included to obtain accurate results. In order to acquire converge result, in the practical 
implementation of QUAPI, one usually need to choose Afc max large enough so that the 
response function reduces to negligible value within the length of Ak max At. However this is 
a hard task, Since augmented propagator tensor A^ kmax "> is a vector of dimension (M 2 ) Afemai 
(M is the system dimension which is two here), and the corresponding tensor propagator 
rp{2&k max ) j s a ma trix of dimension (M 2 ) 2Ahmax , the QUAPI scheme scales exponentially with 
the parameter A/c max . Thus one can not proceed the QUAPI calculation with very large 
A& m ax, and usually Afc max is chosen less than 10 for M = 2, and even smaller for larger M. 



In summary, one have to se 
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result. As discussed in 

of At =0.25/Q (Jl is the tunneling rate of the 
stability for zero bias (as shown in Fig. 2 in 



ect appropriate At and Ak max to achieve stable and accurate 
21| for the SBM problem with the ohmic bath, with the choice 
Dare system there), only Ak max = 5, 7 reaches 



20]) and only Ak max = 7,9 reaches the long 



time limit for nonzero bias (as shown in Fig. 5 in |20|). For our system, it is possible for the 
implementation of QUAPI only when V is large enough, since the memory time is too long 
to implement QUAPI when V is small. 

B. QUAPI2: treatment from the 2nd point of view 

As discussed above, the direct application of QUAPI method is impossible when V is 
small, how ever, we can tackle this problem from the second point of view, that is qubit 
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is first coupled to a harmonic oscillator, which itself coupled to a Ohmic bath. Since the 
memory length of the Ohmic bath extend over only a few slices of the time steps. It enable 
the use of the QUAPI if we take the over all system of qubit-HO as the H and the Ohmic 
bath as the environment. It work well especially for the weak F in this case, since the non- 
adiabatic effect is weak for the small coupling and H provides a reasonable zeroth-order 
approximation to the dynamics. Therefore At can be fairly large, and the result converges 
quickly as Ak ma x increases. 

However, since the qubit-HO system possess an infinite number of energy levels, it prohibit 
the direct use of the QUAPI. We have to truncate the qubit-HO system into a smaller sub- 



space for the practical implementation of the QUAPI method. Similarly to the Ref. [22| . 
we first diagonalize the qubit-HO system in a large N dimensional space and then only 
preserve the lowest M energy eigenstates, consequently the system coordinate operator X is 
also truncated to the M dimension operator Xm in the eigen energy representation. Then 
we diagonalize the Xm to transform into the DVR basis. Since now the system is only M 
dimensional, the implementation of QUAPI method becomes feasible. 

Parameters in QUAPI2: (i) At and (ii) Ak max are the same as QUAPI1. And two 
additional parameter appears in QUAPI2: 

(iii) One parameter is the dimension N of the Hilbert space of the qubit-HO system. The 
Hamiltonian of the qubit-HO system can be numerically diagonalized in the N dimension 
of the Hilbert space. N is kept fixed as N = 400 in our calculation, since it is big enough 
dimension for the diagonalization of the qubit-HO system. 

(iv) Meanwhile, we have employed a second parameter M, which is the lowest energy sub- 
space of the dimension N of the Hilbert space of the qubit-HO system. We first diagonalize 
the qubit-HO space in the larger dimension N of the Hilbert space, to get more accurate 
low energy eigenstates and calculate the physical quantities in the M dimension subspace 
with less numerical effort. Here we should choose larger M for stronger qubit-HO coupling 
9o- 

V. RESULTS AND DISCUSSION 

According to Eq. ©-© and Eq. (fT3]) -(TT5 ]) . P(t) is obtained according to TRWA method. 
Here we report P(t) as a function of time in Fig. 2 and Fig. 3 for the off-resonance case 
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(A = O.lfi), and in Fig. 4 and Fig. 5 for the on-resonance case (A = Q). For the off- resonance 
case as shown in Fig. 2 and Fig. 3, the decoherence is always enhanced by increasing the 
HO-boson coupling T no matter that 7rrA is larger or smaller than g . However, for the 
on-resonance case, the decoherence is enhanced with increasing T when 7rrA > g as shown 
in Fig. 4, Whereas reduced with T when -nTA < g as shown in Fig. 5. 

To check the peculiar results, we also calculate the population difference P(t) by QUAPI 
method. For 7rrA > go we do the QUAPI from the first point of view (QUAPI1) and the 
result is reported in Fig. 3 and Fig. 5, the time splitting is set fixed as At = 0.6/A with 
varying memory steps Afc max . In Fig. 5, one can find that when T is large, i.e. T = 0.4, 0.5, it 
is easy for the calculation to converge and the result is good enough within Afc max = 3. As T 
decreases, it become harder for the evaluation because of the long memory time. In order to 
converge to our analytical result, it needs AA; max = 5 for T = 0.3, AA; max = 7 for T = 0.2, and 
even higher Afc max for T = 0.1 (which almost runs out our numerical resources). For even 
smaller T, such as T = 0.02, the evaluation from this point of view is practically impossible. 
Therefore, for 7rrA < go we do the QUAPI from the second point of view (QUAPI2), which 
is reported in Fig. 2 and Fig. 4. Here the time splitting is set fixed as At = 0.15/ A in Fig. 2 
and At = 0.3/A in Fig. 4, the dimension for diagonalization of TL-HO system is N = 400 
and truncated dimension for QUAPI2 is M = 2 in Fig. 2 and M = 6 in Fig. 4. One can see 
that all these results converges to our TRWA results. 

To understand the qubit behavior, we can explore the damping rate according to the 
TRWA method. One can find that in the near resonance case (A pa Q), level repulsion 
occurs, two characteristic frequencies dominate the qubit dynamics. When the HO-bath 
coupling is weak, the result should agree with that of the Jaynes-Cummings model, that is 
the peak frequency oo p ~ A ± g {(Jl]. Therefore, from Eq. (fTS"]) . we get 

* K S+fc < 21 > 

where we have approximated 77 ~ 1, uj p + r] A w 2A and u p + Q ~ 2A. Therefore, one expect 
7 P oc T for 7iTA<t^go and 7 P oc ^ for nTA^go, which is in accordance with the numerical 
results. Admittedly, as T becomes larger, the above mentioned analysis is not a accurate, 
since the frequency shift of u p will become much more complex because of the dressing of 
phonons. However, from the numerical result, one can find the analysis captures the main 
physics of the coupling dependent behavior. 
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VI. CONCLUSION 



In conclusion, The non-Markovian dynamics of a qubit under the decoherence of struc- 
tured environment is investigated without RWA and Markov approximation. One point of 
view of the problem is the spin-boson model with a Lorentz shaped spectral density. An 
alternative view is a qubit coupled to harmonic oscillator (HO), which in turn coupled to 
a Ohmic environment. Two different methods are applied and compared for this problem. 
One is a TRWA method which is mainly analytical and the other one is the QUAPI method 
which is a numerical scheme based on an exact methodology. The TRWA method can be 
applied from the first point of view. And the QUAPI method can applied from both points 
of views and called QUAPI1 and QUAPI2 respectively. QUAPI1 only works well for large 
T. Since the memory time is too long for the practical evaluation of QUAPI when F is 
small. And QUAPI2 works well for small T, since the non-adiabatic effect become more 
important as F increases, consequently smaller time-step At and more memory steps Ak max 
are required to obtain accurate result which also quickly runs out the computational re- 
sources. We find that the TRWA method works well for the whole parameter range of F 
and show good agreement with QUAPI1 and QUAPI2. On the other hand, we find that the 
decoherence of the qubit can be reduced with increasing coupling between HO and bath, 
which may be relevant to the design of quantum computer. 
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Figures Captions 

Fig. 1: Real and imaginary parts of the bath response function for the Lorentzian spectral 
density. The memory time decreases to a finite range as T increases. 

Fig. 2: The population difference P(t) as a function of time for the off-resonance case 
(A = 0.1O), the parameters are go = 0.1A and 7rrA < g , the QUAPI2 parameters are 
N = 400, M = 2, and At = 0.15/A. 

Fig. 3: The population difference P(t) as a function of time for the off-resonance case 
(A = 0.1Q), the parameters are go = 0.1 A and 7rrA > g , the QUAPI1 parameter At = 
0.6/A. 

Fig. 4: The population difference P(t) as a function of time for the on-resonance case 
(A = fi), the parameters are go = 0.1A and 7rrA < g , the QUAPI2 parameters are 
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N = 400, M = 6, and At = 0.3/A. 

Fig. 5: The population difference P(t) as a function of time for the on-resonance case 
(A = Q), the parameters are g = 0.1A and 71TA > g , the parameters are A = f2, g = 0.1A 
and the QUAPI1 parameter At = 0.6/ A. 
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FIG. 1: Fig. 1 
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FIG. 2: Fig. 2 
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FIG. 3: Fig. 3 
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FIG. 4: Fig. 4 
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